Coupled and Extended Quintessence: theoretical differences and structure formation 



Valeria Pettorino^, Carlo Baccigalupi^ 
^ Institut fiir Theoretische Physik, Universitat Heidelberg, Philosophenweg 16, D-69120 Heidelberg, Germany. 

2 SISSA/ISAS, Via Beirut 4, 1-34014 Trieste, and 
INFN, Sezione di Trieste, Via Valerio 2, 1-34127 Trieste, Italy. 

The case of a coupling between dark energy and matter (Coupled Quintessence) or gravity (Ex- 
tended Quintessence) has recently attracted a deep interest and has been widely investigated both in 
the Einstein and in the Jordan frames (EF, JF), within scalar tensor theories. Focusing on the sim- 
plest models proposed so far, in this paper we study the relation existing between the two scenarios, 
. . . isolating the Weyl scaling which allows to express them in the EF and JF. Moreover, we perform 

OO ' a comparative study of the behavior of linear perturbations in both scenarios, which turn out to 

, behave in a markedly different way. In particular, while the clustering is enhanced in the considered 

■ CQ models with respect to the corresponding Quintessence ones where the coupling is absent and to 
the ordinary cosmologies with a Cosmological Constant and Cold Dark Matter (ACDM), structures 
in EQ models may grow slower. This is likely to have direct consequences on the inner properties 

■ of non-linear structures, like cluster concentration, as well as on the weak lensing shear on large 
scales. Finally, we specialize our study for interfacing linear dynamics and N-body simulations in 
these cosmologies, giving a recipe for the corrections to be included in N-body codes in order to take 
into account the modifications to the expansion rate, growth of structures, and strength of gravity. 
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I. INTRODUCTION 



The increasing amount of cosmological observations has been recently pointing out that the Universe is nearly 
spatially flat, with an expansion rate of about 70 km/s/Mpc and with a 76% contribution to the total energy density 
due to a non-clustered, negative pressure component responsible for the acceleration in the cosmic expansion, known 
^ ' as dark energy. The energy level of this component is markedly lower than almost all known physical quantities, except 
perhaps the neutrino mass differences, and this motivated big efforts in order to understand its physical properties. 
Interesting hints have been obtained studying scenarios in which the dark energy is dynamic, i.e. described by an 
evolving and fluctuating scalar field, the Quintessence. In particular, several investigations have been devoted to 
' study the possibility that the dark energy might be coupled to other entities in the Universe, altering the dynamics 
00 of the Quintessence field itself as well as the evolution of structure formation. It is then essential to understand the 
impact of such a coupling on both cosmological expansion and perturbation dynamics, in view of a clear prediction 
of observable effects for the next to come experiments. 

Two main efforts have been recently carried out in a nearly parallel way. On one side, within the usual framework 
of General Relativity (Einstein Frame, EF) the effect of a coupling between dark energy and dark matter (Coupled 
OO , Quintessence, CQ) has been deeply investigated both in the linear [l| and non- linear regimes 12, J3|, including the 
' effects of such models on the Cosmic Microwave Background (CMB) [3, SQ, on Supernovae 0, HI^Ij on the matter 
power spectrum fiol, o n structure formation [111, , also via N-body simulation [lj| , as well as different choices of 
the coupling ili, E El, [13, E [ill ■ The possibility of an interaction between dark energy and neutrinos has also been 
. , investigated [20 , 121| , leading recently to the proposal of g rowing matter scenarios [12, [H, [l^l . Instabilities within 
. coupled dark energy theories have also been addressed |25l. l26l \2l\ . 
" " On the other side, a wide variety of analysis has been carried out for evaluating the possibility that the gravity 
perceived by dark energy and matter deviates from General Relativity (GR). In this latter case, extensions of GR 
have been considered, in which the dark energy might derive from a non-minimal interaction to gravity via an explicit 
couplingbetween Quintessence and the Ricci scalar (Jordan Frame, JF). This is the case of scalar tensor theories 
[H m,!!!, [m, [H, [M m, HE [H, [13, [H, [H, [13], which are known as Extended Quintessence (EQ) scenarios in the 
framework of dark energy. One of the most interesting results in this framework consists in the so called 'gravitational 
dragging' [36j : the coupling to gravity introduces an effective potential in which the quintessence scalar field rolls, 
acquiring a dynamics which is a purely gravitational effect. Eventually, EQ energy density perturbations can be 
dragged to cluster through the coupling with the gravitational potential activated by the dark matter, its density 
perturbations reaching non-linearity at scales and redshifts relevant for the structure formation process. In particular, 
the growth of dark energy perturbations in the non-linear regime within EQ has been addressed in (40j . 
Furthermore, in the case in which, in the EF, dark energy couples universally to all matter fields, the two parallel 
investigations can also be conformally related by Weyl scaling, which allows to go univoquely from one representation 
to the other. Although the resulting physical effects must be equivalently present in both frames, some features might 
more naturally appear in one representation than in the other and viceversa. 
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In this paper, we identify the theoretical relation between the simplest and most popular examples of both coupling 
choices for CQ and EQ. On the basis of this analysis, we work out the main phenomenological differences in terms of 
the growth of linear perturbations, which turn out to behave in a markedly different way for the two models considered. 
Finally, we propagate our investigation to the inputs required by N-body simulations in view of future applications, 
giving a table of the modifications required in both scenarios for taking into account the modified growth rate of 
perturbations and the strength of gravity. This work is organized as follows: we first devote sec.(jn| to Weyl scaling, 
recalling how it is possible to relate EQ to CQ via a conformal transformation; afterwards, in fsec lllip we investigate 
the background dynamics, comparing CQ in the simple case of a constant coupling (jlll Ap and EQ in scalar tensor 
theories (IIII B[) : we then proceed by studying linear perturbations fsec lIVj) . again presented within both cosmological 
approaches pV Al HVBp . We will then point out in particular, how in the newtonian limit (|V]) the Euler equation is 
modified (sec IV Al sec lVBp . in view of N-body simulations to be realized in both frameworks and for which we specify 
the required corrections fsec IV Cp . Finally in (sec lVip we compare the results and set our conclusions. 



II. WEYL SCALING 



Before handling a parallel analysis of the background and linear perturbations in CQ and EQ models, it is compelling 
to clarify the fact that the two frameworks are strictly related through a conformal transformation called Weyl scaling 
[41I |4^. |43|. [4^ [isl. l46j . It is important to stress that altering General Relativity via scalar tensor theory (Jordan frame, 
JF) is equivalent to coupling a scalar field universally with all matter fields within General Relativity (Einstein frame, 
EE). Although the two reference frames must lead to identical observable effects, the description of the same model 
can in fact rely on different and unequally simple equations when seen in the two frames. In this sense, simplicity 
and naturalness of one frame are not in general preserved in the other one. Going from one frame to the other and 
viceversa according to the model, profiting of this frame equivalence, can therefore represent a valuable tool to easily 
explore a wider set of choices for the coupling functions and the potentials involved: a nasty choice in one frame 
might, in other words, look much simpler in the other frame, allowing us to easily stress remarkably new features, 
more difficult to achieve in the frame of origin. A frame-independent approach has also been considered [45j leading to 
a reformulation of the Boltzmann equation and linear perturbation theory in terms of frame-independent quantities. 
Authors usually refer to CQ models when they consider a quintessence scalar field coupled, in the Einstein frame, to 
the dark matter fields. The action considered in this case usually appears as follows: 



S 



(1) 



where g is the determinant of the background metric, R is the Ricci scalar, k — SttG where G represents the 
gravitational constant, as found in Cavendish like experiments. In the assumption of flat Friedmann Robertson 
Walker (FRW) cosmologies, the line element can be written as ds^ — — dt^ ■\- a'^{t)5ijdx^dx^ , where a{t) is the scale 
factor and t represents the cosmic time variable. The choice of m{(j3) specifies the coupling to ^ matter fields while 
'Ckin,i/' includes kinetic contributions from all components different from (j). 

In the Jordan frame, a scalar tensor theory in which EQ holds is in general described by the following action: 



d x\J —g 



^/(0, R) - iz(0)0'^0;^ - C/(0) - moV^V' + -Ckin.v- 

A\^ A 



(2) 



The scalar field (/> characterizing the generalized gravitational interaction has its kinetic energy and potential specified 
by Z{<^) and V((\)), respectively. In this case the mass mo is a constant and does not depend on the scalar field 
Note that we consider here natural units, c = 1. Compared to general relativity, the Lagrangian has been generalized 
by introducing an explicit coupling between the Ricci scalar and the scalar field, achieved by replacing the usual 
ricci scalar R with the function f{ip,R). This new term, which has the effect of introducing a spacetime dependent 
gravitational constant, may either be interpreted as an explicit coupling between the quintessence field (j) and gravity 
(or equivalently, in the Einstein frame, between dark energy and matter), or as a pure geometrical modification 
of general relativity admitting a non- linear dependence on R [s^. For simplicity, the classes of theories in which 
/((/), i?) assumes the simple form f{(j),R)/2 — KF{(j))R/2 have been widely considered in the context of dark energy 
cosmologies. 

Weyl scaling consists in a conformal transformation of the metric which, joined to a redefinition of matter fields, 
allows to rewrite action ([2]) into ([T]) or viceversa: 

ffM'' = A^{4>)9tJ.i^ (3) 
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R = [A(0)]-2{i?_6g^''(lnA);,(lnA);J 



(4) 
(5) 
(6) 
(7) 



where we have used the ~ to identify quantities in the EF and distinguish them from those in the JF. Note also that 
the scaling factor ^(0) is related to the coupling via the following relation: 



1 



(8) 



When applying Weyl scaling to the action ([2|) with Z{(f))\ for simplicity, we obtain that the rescaled action appears 
to be 



S — d X\J~g 



kin,^ 



where 



Z((/.)=A^(0) 1 + - 



6 ^2 



(9) 



(10) 



^(0)=A4(0)C/(0) , 
and both the mass and the kinetic term for matter fields now also depend on 

m{(f>) = moA{(j>) 



(11) 



(12) 



2A(0) 



(13) 



Note that the non standard kinetic term Z in (jlO[) is always positive and one could further re-absorb it by redefining 
the scalar field (j) implicitly into Lp in such a way that: 



I? ) - m 



This allows to rewrite action (151) in the usual form: 



S = / (Tx^/^ 



1 ^ 
2k 



(14) 



(15) 



We have now all the means to find out, given a certain action ([2|) in the Jordan frame, which is the equivalent action 
(|15p in the Einstein frame. In the parallel analysis that will follow in this paper, we will eventually refer to two 
particular choices of the coupling in CQ and EQ scenarios. In particular, we will consider two of the simplest and 
most popular choices, that is to say a quadratic coupling in the Jordan frame for EQ [s^, 113, S, and an exponential 
coupling in the Einstein frame for CQ [T, '2']. In order to have a feeling of a possible comparison between the two 
different models, we illustrate here how they look like in the same, Einstein, frame, after Weyl scaling a quadratic EQ 
model. 

In the case of EQ with a quadratic coupling, we will work with so called 'non-minimally coupled' theories, widely used 
in previous works on this topic (see 34,] and references therein), in which F((j)) is the sum of a dominant constant 
term plus a piece depending on (jj: 



8nG* 



(16) 
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Here k* = SttG*, where G* represents the 'bare' gravitational constant [42|, which is in general different from 
G and is set in such a way that locally 1/k + £,<l>iocai ~ l/(87rG) in order to match local constraints on General 
Relativity. determined by a combination of G* and the local value of f{(f),R) such that G = G,,y \iocai- order to 
recover the equivalent action in the Einstein frame, we apply equations ^ [TOl El fT^ . thus obtaining action ^ in 



which: 



m{c^) = , . (19) 

G 



Notably, in this frame the coupling to matter fields, the kinetic term and the potential do not look as straightforward 
as in the original action. 

Another very popular choice is represented by Induced Gravity (IG) [H, [H, |43, Sll, in which only the quadratic 
coupling 

F(^) = (20) 

is considered and no constant term is present. It can be easily seen, by applying equations ([51 [TUl [TTl [T^ that in this 
case we recover action ([9]) in which: 

1 



^(<^) = 3^(1 + 60 , (21) 



m(0) = . (23) 

In this case it is also straightforward to solve equation (fH)) analytically, obtaining 

= De'^''^ , (24) 

for > 0, where G and _D are constants and we have defined G = ^/^/{l + QO- The latter expression, substituted in 
(l22l [23| allows to recover action (fT5|) in which: 

C/(^) = ^e-^'^U{H^)) , (25) 

m(^) = ^e-^-. (26) 

Applying Weyl scaling to IG theories, leads, as expected [H, H^, to CQ models in which the mass of matter fields 
depends exponentially on the scalar field ip. We therefore recover the simplest and most popular choice in CQ 
investigations flt| and we will ourself adopt this choice in this paper when dealing with CQ in the Einstein frame. 
Furthermore, since U{(j)) is usually almost constant in order to match the present constraints on the amount of 
Quintessence, the transformed potential in the Einstein frame behaves as an exponential. 

As appears clear from the previous analysis, IG theories in the Jordan frame (or equivalently CQ models in the Einstein 
frame with an exponential coupling in the Lagrangian) have a crucial difference from the non-minimal coupling choice 
described by eq. pB]) in the Jordan frame. Indeed, IG (or CQ with an exponential coupling) forces the coupling 
constant f to be positive in order to get the right sign for attractive gravity in action Q ; on the contrary, the choice 
(fT6)) also allows negative values of ^ and can therefore lead to an opposite correction to background and perturbation 
dynamics, as we will see in detail in the following section. 



III. BACKGROUND DYNAMICS 



In this Section we define the CQ and EQ models we study in the following, and review their dynamics for what 
concerns the cosmic expansion, highlighting the aspects which will be relevant for the behavior of linear perturbations. 
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A. Coupled dark energy 



We will follow as much as possible previous works on this topic [l|, |3| , considering the simplest models in which the 
dark energy scalar field and CDM are coupled via a constant coupling appearing in the conservation equations. 
The equations used to describe the background evolution of each component (a) involved in the interaction follow 
from the consideration [i^ that the coupling can be treated as an external source acting on each stress energy tensor 
T^u{a) ill such a way that the total stress energy tensor is conserved: 

Ti:,, (e) = (27) 



(0) = C-.T.^^, , (28) 

where Cc is a constant, Tc is the trace of T'^^(^c) and the subscripts (c), (0) stand for CDM and dark energy respectively. 
Baryons do not couple with dark matter or dark energy. The constant coupling term here used can be achieved starting 
from the following Lagrangian: 



1 „ 1 
2k 



(29) 



in which the mass of cold dark matter fields depends exponentially on (j) |58j : 

to(0) = moe^^"'^ , (30) 

thus corresponding, in the JF, to IG cosmologies, as seen at the end of the previous section. Here, for simplicity, we 
do not include the baryonic component. 

Equations (j27l28p specialized for v = 0, provide the conservation equations for the energy densities of each component, 
given by p'^ — —3Hpc — CcPc4>' and p'^ = —ZTLh^ + CcPc4>' for the two coupled species and = —^Hpr since Cr = 
for radiation, pr takes into account the radiation contribution, including both photons and neutrinos, both decoupled 
from the other species. We have further defined h^f, = Pa + Pa = (j) ^/o^ where the latter equality comes from the 
expression of the energy density p^ defined from the (0, 0) component of the stress energy tensor present in (|28p and 
therefore given by the usual 

P0 = ^^ + C/(0), (31) 

where U {(f) is the potential in which the dark energy scalar field rolls. Note that is still equal to the standard one 
characterizing uncoupled quintessence models since, as we have seen, the coupling is provided in (|27[|28p as an external 
contribution to the standard stress energy tensor; nevertheless, unlike what happens in the uncoupled quintessence 
case, only the sum pc + is conserved while p^ is not conserved by itself. 
Analogously, the pressure of the quintessence scalar field is defined as: 



For what concerns the background dynamics, the evolution of the scalar field is described by the Klein Gordon 
equation, which reads as follows: 

0" + 2n4)' + = a^CcPc , (33) 

where U^^ is the derivative of the potential with respect to 4>. 
The CDM conservation equation can be formally integrated out 

^ P^f.-cA4-<t>o) (34) 



a 



giving the dependence of pc from the scale parameter, once that the dynamics of 4>{a) is known. 
The evolution in time of the scale factor is provided by the usual Friedmann equation, so that: 

■H^^a' P^ + P^ + P^ . (35) 
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For convenience we have defined p = SirGp = p/M"^ and (p = (^jM in units of the reduced Planck mass = I/SttG = 
Mp^ /Sit where G is the gravitational constant. The coupling constant C = CM is therefore dimensionless; we have 
then omitted the " and redefined p = p, (j) = (j) and C = C for simplicity. 

A detailed analysis of the background has already been addressed in [l| in the case in which U{(f>) is an exponential 
and in [1, for a more general treatment. In particular, the critical points of the trajectories in the phase space 

of the system have been widely investigated in We just recall here that in the case in which both the potential 
U{(l)) and radiation are neglected, it is straightforward to find the critical points of the system, corresponding to 
d(j)crit/da = (— -s/H, +VQ, '^Cc) where we have defined a = logja(r)]. The third critical point is the only one 
providing a stable attractor during matter dominated era (MDE) [l| ^ along which the evolution of the dark energy 
scalar field follows the trajectory 



</>(a) = 2Cc log + 00 , (36) 
where we have put todays value of the scale parameter equal to oq = 1. Substituting in ([M)) . we get: 

Pc = -g-a " ■ (37) 

We can now make use of the Friedmann equation, specified in the case in which there is no radiation and the potential 
is negligible, to obtain the scale parameter dependence on the conformal time during MDE: 



aMDEir) OCT I - "i (38) 

as well as 7i(r): 

2 1 

UMDEir) = Y^2C^T ' ^^^^ 

or equivalently, in terms of the cosmic time aMDE{t) oc t[2/(3+2C^)] HMDsit) = 2/[(3 + 2C^)t]. We stress here 
that, as we can see from the four latter expressions, the effect of the coupling on the Hubble parameter and on the 
scale factor during MDE, when the potential is negligible, does not depend on the sign of the constant Cc, both of 
them being only a function of C^. With this choice of the coupling, we can thus only have a one direction effect, 
being that of slowing down the universe expansion with respect to the standard case in which Cc — 0. 
We can numerically solve the equations describing the background, in order to obtain the full evolution of the energy 
densities of both cold dark matter and the scalar field, as shown in fig.([T]). In the left panel, the energy densities of a 
cosmological constant (dotted), radiation (long dashed), CDM (solid) and coupled quintessence (dashed) in the case 
in which C = 0.05 have been plotted. We have used an inverse power law shape for the potential Uicj)) which was 
first proposed in [s^l and often used in quintessence models since then [5l|, [s^l : 

^0 



Ui^) = Uo[jjj , (40) 

where (f>o is the value of the scalar field today and A will be typically fixed to ~ 0.2 in order to have a flat potential and 
a consequent equation of state for dark energy which is close to —1 today, in order to be consistent with observations 
[53I [5^ . Issj . Another popular choice is given by an exponential potential [l|, [s^, HI] , though we are not interested 
here in the explicit expression of U{4>), which does not modify our considerations. The amplitude of the potential 
()40|) and the initial amount of CDM are fixed in order to have their values today again in agreement with the present 
constraints [H, |H, HI]. As shown in fig.© (left panel), the presence of the coupling has a drastic effect on the 
behavior of the scalar field energy density. In the uncoupled case quintessence behaves pretty much as a cosmological 
constant due to the rather flat potential. If however the coupling to CDM is active, accelerates and the dark energy 
density is attracted by the evolution of the species to which it is coupled: in particular, after a period of transition 
due to the particular choices of the initial conditions, is dragged by CDM to follow a trajectory which strictly 
mimics the CDM path, therefore behaving as an effective pressureless component. We will refer to this effect as 'dark 
matter dragging' in analogy to the 'gravitational dragging' flrst discussed in [3^ in which the coupling of the scalar 
field to gravity via a modification of General Relativity imprinted a similar pattern to the scalar field dynamics, both 



^ The three critical points correspond to the last three critical points in table 1 in [J] . 




FIG. 1: Energy densities of CDM and dark energy. The left panel shows the case of coupled dark energy for C — 0.05 as well 
as for C = vs 1 + z, where z is the redshift. The right panel shows the case of extended quintessence for both positive and 
negative coupling corresponding to ljjbdo ~ 30. During MDE the energy density of the scalar field (the non conserved one) is 
more enhanced in the case of coupled dark energy than in extended quintessence, in which the dynamics of the R-boost has its 
major influence during RDE. The plotted densities are in units of Mpc~^ . 



on background and on perturbations. In this scenario, the matter scaling is also modified by the coupling with the 
field; the resulting regime of cosmological expansion driven by the field and CDM tracking each other, is known as 
0CDM era, see j3| and references therein. 

In both models here considered, (f> rolls down an effective potential provided by the coupling to dark matter (right 
hand side of eql33|) or by modified gravity (right hand side of equation |44l next paragraph). Note however that 
despite the analogy with gravitational dragging, during dark matter dragging the dark energy scalar field does not 
exactly slow roll on the effective potential, as it happens in the gravitational dragging case: when the coupling to 
dark matter is considered, (j) acquires a strong acceleration from the dark matter field, that makes terms like (j)" give 
a non negligible contribution in equation psp, an occurrence which is absent in the EQ case, as we now discuss. 



B. Extended quintessence 

As illustrated in the first section, scalar tensor theories are described by action ^ that we rewrite here for conve- 
nience: 



S 



-!-/(0, R) - ~U{(I))- moipip + £kin,v 



(41) 



We have chosen Z{(f)) = 1 to simplify the notation. We further set f{4>,R) = KF{(j>)R so that the function F{(j>) fully 
specifies the explicit coupling to gravity via its dependence on some scalar field (j>. When the additional scalar field 
introduced in the coupling plays today the role of the dark energy field, leading to an accelerated expansion of the 
universe, we refer to it as extended quintessence, see j37| and references therein. The limit of general relativity is 
achieved with ljjbd ^ 1 where lujbd is defined as 

wjBD = ' (42) 

and is the derivative of F with respect to (p. The current constraints set the present value lujbdo to be larger than 
40000 in the Solar System [63|, and greater than about 120 from cosmology [6l|; being on markedly different spacetime 
scales, the two are considered as independent and complementary [b^]. Within these theories, the Friedmann equation 
describing the expansion of the Universe is modified with respect to the usual case: 



a 



Id)' SHF' 
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Besides the extra term at the end of p3|) . which is usuaUy quantitatively neghgible, the most important effect to point 
out is the presence of the muhiplying factor 1/ F substituting the usual SttG; as a consequence, extended quintessence 
behave as theories in which gravity depends on a varying function and not on a gravitational constant anymore. Note 
that now all fluids but (p satisfy the usual conservation equation p' ^ — STiha including CDM, which now scales as 
Pc = Pcq/o!^ ■ The evolution of the scalar field is explicited by the Klein Gordon equation, which again gets contribution 
from the coupling to gravity, enhancing the background dynamics of (j) and giving rise to the R-boost effect illustrated 



where and are the derivatives of the coupling F{(t)) and of the quintessence potential U{(f)) with respect to 
(f). The term on the right hand side has the effect of altering the potential into an effective potential in which 

the scalar field rolls. This effect, responsible for the gravitational dragging ,36j, is analogue to the one produced by 
the right hand side of equation ([33]) due to the coupling to matter fields, within the Einstein frame, as shown in 
fig. H] right panel): here we plot the energy densities of radiation (long dashed), CDM (solid) and EQ in the case 
in which the coupling to gravity is positive (dashed) or negative (dotted). The energy density of the scalar field 
has a similar behavior during RDE (R-boost effect), independently of the sign of ^, while the two patterns detach 
a bit during MDE. The plot has been obtained solving numerically the equations for the background in the case of 
non-minimal coupling, with F{(j)) given by (|16p . Again, the effect of the coupling on the background evolution is 
that of enhancing the amount of dark energy in the past, due to the gravitational dragging set on by the additional 
term in the Klein Gordon equation (j44[) . Though in the quadratic case the gravitational dragging shows in particular 
during RDE, note that the time at which the dragging ends depends on the choice of the model and on the values 
of the coupling parameter. An exponential coupling in the lagrangian would have a bigger effect, lasting up to MDE 
and more recent times [s^. Although the phenomenology and the energy density scaling is clearly analogous for the 
CQ and EQ models plotted in fig. HI), it is important to stress here that the sign of the coupling leads however to a 
different correction in the Hubble expansion parameter. Indeed, while in coupled dark energy the Hubble parameter 
does not depend on the sign of the coupling constant during MDE feq [5^ . in extended quintessence a positive and a 
negative coupling lead to different sign corrections, either enhancing or reducing the Hubble parameter with respect 
to the standard LCDM case. The effect is shown in fig. ([2]): here we plot the bubble parameter versus redshift for 
LCDM (solid), coupled quintessence with C = 0.1 (dotted) and extended quintessence with a positive (long dashed) 
or negative (dashed) coupling. For the chosen coupled quintessence model, with a constant coupling, the bubble 
parameter is bigger in the past than the usual LCDM case, independently of the sign of the coupling constant. The 
chosen extended quintessence case, instead, can lead to both higher and lower values in the past, depending on the 
sign of the coupling. Note also that the switch in sign does not lead to perfectly opposite contributions: when the 
coupling is negative the effect is bigger than in the case of a positive coupling with the same absolute value. A reason 
for this to happen might be that in the first case the extra term in the Klein Gordon equation adds to the usual 
potential, favoring an easier enhancement of the dynamics of the field; on the contrary, for a positive constant, the 
extra term contrasts the effect of 1^(0), making it more difficult to enhance the dynamics of the field. 

Finally, we recall that in this context, the Ricci scalar can be written in terms of the cosmological content of the 
Universe 



where p fluid and p fluid are the energy density and pressure summed up over all possible cosmological components 
but (p. 



Here we perform a comparative analysis of the behavior of linear perturbations in the two scenarios. As we show 
in a moment, the main difference highlighted in the previous section persist at this level, as in the Newtonian limit 
which is outlined below. 



in [H 



(f," + 2n(t>' + a%^ = —F^^R, 



(44) 




(45) 



IV. LINEAR PERTURBATIONS 
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FIG. 2: Hubble parameter vs redshift for CQ and EQ with a non minimal coupling (for positive and negative couplings). The 
case of ACDM is also shown for reference. The value of ^ corresponds to lujbdo ^ 30. H is in units of Mpc~^ . 



A. Coupled dark energy 

The linearized perturbation equations can be derived, following the notation of [49| . from 



(46) 



with = (— aQajO) being the source term for the a— component. In particular, comparing with eq. (|27p and 

eg. ([28)1 ■ for the case of a constant coupling here considered the source is defined as: 



CcPc4>' 



(47) 
(48) 



For V ~ Q one gets the perturbed energy conservation equations, while for v = i (spatial index) we can recover the 
first order Euler equations. In the Newtonian gauge in which the non-diagonal metric perturbations are fixed to zero, 
the linearized conservation equations for the matter component and for the scalar field, respectively, read as follows: 



fc(* + CM) , 



Sp'^ + SHSpc + Pckvc + 3pc*' = 

6p'^ + 3n{Sp4, + Sp^) + kh^v^ + 3h4,^' = Cc{pM' + <!>' 5pc) , 
h^v'^ + (h'^ + ATih4,)v^ = k5p^ - kh^^ + Cckp, 



Sp'^ + AHSpr + k-prVr + 4p,.*' = , 
o 



(49) 
(50) 
(51) 
(52) 

(53) 

(54) 



where k is the wavenumber, Va is the velocity perturbation (along the direction of k) for the a species and $ is 
the gravitational potential coming from the metric perturbation; in particular $ = — if we neglect the anisotropic 
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stress. Note that the latter is rigorously zero for the present scenario, as opposed to the EQ case as we will see in a 
moment. Thus the only source of that here is represented by neutrinos, due to their viscosity for having decoupled 
in a markedly relativistic regime. Also note that the two latter equations hold only if we neglect baryons, otherwise 
Thomson scattering with photons should be taken into account. A more convenient way to express equation (j49p is 
to rewrite it in terms of Sc = Spc/pc as 



S' + kvc + 3*' 



(55) 



The gravitational potential $ is related to the background and perturbed densities and velocities of the various 
components via the following expression: 



-r ^ \ 



a 

2fc2 



Spc + Sp^ + Spr+S"^ [ PcVc + ^PrVr + h^V^ 



Or, equivalently, in a more compact form and in terms of the f^(a) as 

3, 



^a,(<5, + 3(l+c^,0Aw, 



(56) 



(57) 



where we have defined A = Ti/k and where Sic = {o? Pc) / {i'H'^) and the sum is extended to all components, including 
dark energy. Note that dark energy perturbations contribute to the gravitational potential and therefore on the dark 
matter perturbations through eg. ([55)1 . 

Deriving with respect to the conformal time, making use of equations ([^ [5D1 [STl [5^ as well as of the background 
conservation equations and substituting the expression of $ given by (|56p we get the following simplified expression 
for $': 



2k 



The energy density perturbation for the dark energy scalar field can be expressed in terms of 
derivatives in the following way: 

5p^ = p^8^ = 4 yiH)' + a^U,^5c^ + ((/.')'*] ■ 



(58) 
and their 

(59) 



We recall that the latter expression, holding in the newtonian gauge, can be derived from (5Tq(0-) of the dark energy 
scalar field when comparing it to the usual one for a fluid of pressure and energy density p^ (or, equivalently, 
directly perturbing pip V Analogously, for the velocity perturbation in the newtonian gauge, one has 



(60) 



which can be derived from the ST'^f^^^ component of the perturbed stress energy tensor. Making use of the two latter 
expressions, $ and can be rewritten as 



2fc2 - (0') 



l\2 



(61) 



= - 



27i* + 3- 



k 



(62) 



Note also that due to the definition of the energy density and pressure of the quintessence scalar field (|31] [32l) the 
perturbed pressure appearing in eg. ([5^ is simply equal to: 



2Ua 



(63) 



As regard to the dark energy scalar field perturbation, 50 can be either written using (|59p or it can be gained when 
solving the perturbed Klein Gordon equation: 



b" + 2H(50' + (fc2 + a^U^^^)6<j) + 40'*' - 2a^U,^^ = Sn^ileCe [6c - 2*] 



(64) 
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FIG. 3: Growth factor as a function of redshift for coupled dark energy and extended quintessence with a non minimal coupling 
(for positive and negative couplings). The case of LCDM is also shown for reference. The value of ^ corresponds to lojbdo ''^ 30. 



In fig.® the density perturbation of CDM devided by the scale parameter and normalized to one today is shown; the 
latter quantity is the cosmological growth factor for CDM. The curves have been obtained by numerical integration 
and correspond to the case of CQ (dotted) as well as EQ case, for a positive (long dashed) or negative (dashed) 
coupling, as we will expose in detail in the following. The ACDM case is also shown for reference (solid): as it is well 
known, in the latter scenario the growth factor is almost constant during the Matter Dominated Era (MDE, (5c oc a), 
while it increases going back in time due to the presence of radiation. The cosmological parameters of these scenarios 
are as defined in the previous section. We postpone further comments to this figure to the end of this section, when 
we expose the EQ case, too. 

B. Extended quintessence 

CDM variables follow the usual uncoupled equations, written here in newtonian gauge: 

d'^ = -kve - 3*' , (65) 

v'^ = -Hvc + /c* . (66) 

Note that, in contrast with the CQ case, as we will see later, in EQ the two gravitational potentials $ and 4^ are not 
equal in module, since the anisotropic stress is in general different from zero. As we show now, however, this is not 
the major source of difference between the two scenarios. 

As shown in (ssl . [36| the perturbed Einstein equations can still be written formally as 

SGI = k6TI; , (67) 
where is the usual Einstein tensor and STjf is the total stress energy tensor, summing up all species: 

STi; = STI^^f^ + 5r,^(^) . (68) 
In particular the conserved stress energy tensor for cj) is given by 

= 5Tr'm + ^T^^-l^] + STr^'^id^] , (69) 

where, in the newtonian gauge: 

5TI^ "^[0] = ^ - 0'5</)') - U,^54> (70) 



12 



3n 



2* 



5! 



F' + ^JlsF' + 

2 



R 



SF 



jT,^f™"[0] = ^5F- 



(71) 
(72) 



Carrying the last term of expression (|72p to the left hand side in the perturbed Einstein equations, we can rewrite 
eas (|67)) as 

(73) 



FSG^, = STr^i^] + STr-'-'m + ^^SF + <5r,^(/,„,) 
The (0,0) component of ([75)1 reads 



— ^SP =^ — 



^ (HF + F')^-% {2nF + F') 



(74) 



i? 



F.^R 



mSF' - [k^ + iTT - - ] 5F ] + pfSf 



The (0, j) component of ([75)1 gives the following equations: 



^ J jiuid 



(75) 



Combining the two previous equations to get rid of we get the Poisson equation in EQ cosmologies, which is 

/2\ 



k'+6n'- — + 





1 


— * 












5 (n + 


-) 


2Fy 


m F'" 


) 5F 









1 

+2— 



(76) 



3 F' 
^ 2 F 

F.d>R 



2 



where we have used that in the Newtonian gauge STf = vrhr. Finally, the Klein Gordon equation which describes 
the evolution of the EQ scalar field is 

^2 



(50" + 2nS<j)' 



(77) 



2 2 

(*' - 3*') + y (F^i? - 4[/,^) * + ^F^SR , 



where 



(37^* - 3*')' + (37^* - 3*') - (fc^ - m' + 3H^) * - 2k^^ 



(78) 



These equations hold in general for any choice of the function F{(j)). In order to solve the equations numerically, we 
again set the case of non- minimal coupling, in which ([IB)l is valid. 



The results of the numerical integration are shown in fig. ([5]). As a reference, we normalize the growthfactor to a 
common value at present time. In the CQ case, the coupling to CDM has the effect of lowering the growthfactor 
with respect to the ACDM case; this means that for a fixed primordial normalization of the perturbations, the CQ 
structure formation is enhanced with respect to a ACDM case, independently on the sign of the coupling constant, 
as it is known from earlier works The EQ phenomenology is completely different, as the structure formation 
may be slower than the ACDM case, depending on the sign and magnitude of the coupling constant. This may have 
interesting consequences for constraining these theories from the observations of cosmological structures in the mildly 
and full non-linear regimes, as we see in the next section. Note that, consistently with the analysis of the previous 
section, in the EQ case the departure from the ACDM case is again bigger for a negative coupling. 
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V. NEWTONIAN LIMIT 



A coupling between dark energy and matter or gravity could indeed have observable effects on structure forma- 
tion, due either to the interaction between species which behave differently with respect to gravity or because of a 
modification of gravity itself. We perform here a derivation of the Newtonian limit for cosmological perturbations in 
CQ and EQ scenarios [1, [s^, |3^ . This part of the work represents the interface to numerical simulations of structure 
formation in these scenarios. Previous works on this subject considered the cases of CQ only [3|, 113|. We do expect 
to be able to track the differences outlined in the previous sections also in the present context. 



A. Coupled dark energy 

We strictly follow the notation of [1, [l^ and we specialize the previous general equations to the case in which 
A = Ti/k ^ 1. This choice corresponds to the Newtonian limit, that is to say scales much smaller than the horizon. 
When the newtonian limit holds, the expressions of both the gravitational potential 4" and its derivative can be 
largely simplified and become approximately equal to: 

* = -L [^'{5(j})' + m^'5(j} + a^u,^S(j) + s^cn^Sc] , (79) 
^' = -^(f>'S(i)-n^ . (80) 

Substituting in the perturbed Klein Gordon equation ([64]). we get the following equation describing the evolution 
of <j): 

5(j)" + 2^(50' + (fc^ + a2f/,00 - 2(j)''^)5(j) + *(-47i:0' - 2a'^U,^ + QTi^VLcCc) = BH^r^cCc^c • (81) 

In the Newtonian limit k ^ Ti. and the latter equation can be further simplified: the term containing # can be 
neglected since $ is order fc^^; within the terms multiplying Sip , we neglect the contribution of the potential, which 
has effect only at very recent times, when the acceleration takes over; furthermore, the remaining piece k^ — 2(f>''^ can 
be rewritten as fc^(l — 2X'^{d(j)/da)'^) ~ fc^ where a = log a. Supposing that also Scf)" and Sip' can be neglected (this 
can be further checked a posteriori), eq. (|8ip in the newtonian limit reduces to: 

S^ - SX^flaCcSc , (82) 

stating that Sep is of order and can therefore be neglected, together with its derivative, in eq. ([79|) with respect to 
the last term: 

* ^ h^n.Sc . (83) 

By looking at the right hand side of eg. ([50]) we can formally define the quantity: 

*c = * + CcS^ , (84) 

which reads 

*e = ^X'nM^ + 2C,') , (85) 

where we have used the approximated expression (|82p for Sip obtained by solving the KG equation in the newtonian 
limit. The latter equation reads in real space as 

V2*e = -ypA(l + 2C2) , (86) 

and if we explicit the units contained in p ^ we then have that $c behaves like a new gravitational potential satisfying: 

V2*e = -inCa^pJc , (87) 



^ Recall that our energy densities are equal to p = SirGp and we have redefined p = p for simplicity. 
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where we have defined G* as 

= G(l + 2Cl) . (88) 

In this way, we formally recover the usual Poisson equation, where the effect of the coupling is partially included in 
the redefinition of the gravitational potential $c and partially in a varying gravitational constant G* . 

The equations for the CDM component (|55l and [50]) can also be analogously simplified and read, in the Newtonian 
limit, as: 

5c = -kvc , (89) 



v', + {n~ Cc4>')vc = -k^o ■ (90) 
Deriving the first equation and combining it with the second, it is straightforward to obtain 

5'^ + {H- Cccj)')s'^ - ^n^nM^ + ^c^c) = o . (91) 

Furthermore, equation (j90p allows us to clarify how the Euler equation, relevant for the interaction of dark matter 
particles within structure formation, changes due to the presence of a coupling between dark matter and dark energy 
(Einstein frame). Taking into account expression (I84p and rewriting (|90p in real space coordinates, the equation finally 
reads as a modified Euler equation: 

Vi;^ + {n- Gc0')Vwc + ^n^ncScii + 2C^) = . (92) 

Supposing now that CDM is concentrated in one particle of mass rric at a distance r from a particle of mass Mc at 
the origin we can rewrite the cold dark matter density contribution as 

where we have used the fact that a non-relativistic particle at position r has a density given by mcnS{r) (where S{r) 
stands for the Dirac distribution) with number density n — e^'^^i't'^'t'o) given by equation ([34)) and we have assumed 
that the density of the Mc mass particle is much larger than pc- The Euler equation in cosmic time dt — adr then 
reads: 

\/vc = -HVvc 5^-^ , 94 

where we included the mass scaling in a new redefinition of the gravitational constant 

G = G*e"^<=("*-"*") = G(l + 2G2)e-^'^('^-'^°) , (95) 
and the further effect of the coupling on the friction term is a redefinition of the Hubble constant: 

H^H^l-Cc^^ . (96) 

Though the choice of including the exponential correction in a redefinition of the gravitational constant might be 
convenient for numerical purposes, we can as well rewrite eq. (l94p as: 

47rG.A?c(5(0) 

Vvc = -HVvc . , 97 

in which gravity is still governed by a gravitational (though modified) constant and the time dependent extra con- 
tribution is included in the mass Mc = Mce"'^^'^'^"'^''-' of dark matter: when dark matter particles interact with each 
other, each dark matter particle sees the other one as possessing a 'modified' mass which takes into account the 
non-standard scenario due to the coupling. 

Along the attractor ((36)) . the full correction to the gravitational constant and to the expansion assume the following 
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d ■'■■'4 



: CQ with C = 0.05 

1-12 - CQ with C = 0,1 

EQ with f = 0.21 



1.1 
1.08 
1.06 
1.04 
1.02 
1 

0.98 
0.96 



EQ with f = -0.19 
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FIG. 4: Gravitational correction as a function of redshift for CQ and EQ models. The value of ^ corresponds to lujedq ~ 30. 



expressions: 



G = G(l + 2C^)a 



(98) 



H = H{l-2Cl) , 



(99) 



(100) 



This clearly shows how the corrections behave regardless of the sign of the coupling constant. At the end of this 
Section, we'll focus on this aspect again, comparing the corrections to the strength of gravity in CQ and EQ. 



B. Extended quintessence 



In the limit of small scales, in which fc 3> 7i we can proceed in simplifying the equations for EQ in a way which is 
analogous to the one adopted in CQ. Neglecting time derivatives (negligible in the newtonian limit with respect to 
spatial derivatives) and the potential, the Poisson equation ([7^ reduces to: 



which reads, in real space, as: 



F 



<^ V72 ; 



(101) 



(102) 



From the traceless part of the perturbed Einstein equations (component (i,j)) we get the gravitational potential 
relation 



(* + *) = 



SF 



(103) 



so that in general the right hand side, representing the anisotropic stress, is different from zero in scalar tensor theories. 
Equation (jl02p can then be equivalently written, in terms of '3/, as 



F 



0V72; 



(104) 
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In the approximation we are considering, the perturbed Klein Gordon equation given by (|77p only receives contribution 
from 

(50 = F0[* + 2*] . (105) 

Substituting ()103p . we obtain: 

FFa, FFa 
5<j)^ -% ^ = i-^^*. (106) 

Further substituting in (jlOip . the Poisson equation can be rewritten in the usual way 

~^ F 

or, in real space, as: 



— *E = -^Pf6f , (107) 



V2*E = -^«'p/^/ , (108) 



if we redefine the gravitational potential as 



Note that we could have equivalently substitute 6(p{^) in (|104p thus obtaining 



F 

if we redefine the gravitational potential as: 



V^*E = ^a'pfSf , (110) 



*E= (l- J— I * . (Ill) 



2F + 2F,^\ 

The Euler equation (j66p is then modified according to: 

v'^ + nva + k\l + — ^^^1*^0. (112) 
V P + P- ^ I 

Substituting the expression for $($e) given by (|109p and making use of the Poisson equation (|107p we get 

2(j^ + 2f|) 
2k'^'''' Fi2F + 3Fl) 

which in real space reads as: 



< + + ^.Pf^f^TW^rvT^TX = , (113) 



3 2(F + 2F^) 

v^^ + nvv, + -n'nfdf ^^^^ ^ J,^^ ^ o , (im) 

where we have defined f^/ = Pf/pcrit and Pcrit = STi-^F/a^. For comparison with the CQ case investigated in the 
previous paragraph, suppose now that the only other fluid present is CDM, concentrated in one particle of mass 
at a distance r from a particle of mass at the origin, such that 

^^M^ (115) 

due to the fact that in this framework cold dark matter is uncoupled from dark energy and approximately scales in 
the usual way pc oc a^^ . Substituting the latter expression in (jll4p we get 
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in terms of the cosmic time, where we have included the extra contribution in the redefinition of the gravitational 
constant, which is now varying in time: 



2{F + 2Fl) 1 
{2F + 3i^2 ) ^^p 



(117) 



The latter formalism is general for any choice of F{(j)). For illustration we specify our final results to the non minimal 
coupling choice given by pB]) . In this case 



G 



1 

SttG. 


Ki- 




1 


1 

SttG, 


Ki- 









(118) 



For small values of the coupling, that is to say ^ <C 1 the latter expression becomes: 

- 1 - SttG^.^" (119) 

which manifestely depends on the sign of the coupling ^. In fig.(l4|) we compare the behavior of the correction to 
the gravitational constant in CQ and EQ theories, for different values and sign of the coupling constants. As we 
stressed in this Section, in the CQ case the correction is independent on the sign of the coupling constant, just as 
expected since the theory corresponds to the case of induced gravity, which, as discussed in the first section, forces 
the sign dependence in order to mantain gravity as an attractive force; note also how also at present the value of the 
gravitational coupling is larger than the corresponding ACDM case, as it is evident from (|98p . On the other hand, in 
the EQ case, the sign has the effect of increasing or decreasing the gravitational strength. 



C. Recipe for Nbody users 



We conclude this section by summarizing the quantities needed as an input for N-body simulations when considering 
CQ models or EQ models as discussed in this paper. An N-body code which wants to take into account the modified 
gravitational strength between different species, as well as the modified perturbation growth rate, needs to read tables 
with the modification listed in table |T] for both sets of cosmologies. The expansion history, represented by the Hubble 
expansion rate, needs also to be modified, although it is not listed in the table. 

Three typologies of corrections are considered: the gravitational interaction between dark matter (DM) particles, the 
friction term due to a modified expansion and the mass of the interacting species. 

When dealing with CQ models, all three corrections are active: the gravitational interaction between DM-DM particles 
changes, the friction term provides an extra contribution due to the coupling and the mass of the DM particles also 
changes with time and depending of the coupling. In the case in which we also consider baryons (B), non coupled to 
the scalar field (j) or to CDM, the gravitational interaction between DM-B or B-B is the usual Newtonian one, since 
non coupled baryons satisfy the standard Euler equation. Furthermore, no correction to the baryon mass is required. 
On the contrary, when EQ cosmologies are considered, only the gravitational correction needs to be applied. The 
effect of the coupling is all included in an effective gravitational 'constant', which changes with time according to 
eq. ljllSp for all the species considered. Both the friction term and the mass term remain unaltered. 



VI. CONCLUSIONS 



In this paper we performed a comparative analysis of the theoretical framework and phenomenology of cosmological 
perturbations in two of the main scenarios where the dark energy is explicitely coupled with dark matter (Coupled 
Quintessence, CQ) or gravity (Extended Quintessence, EQ). Weyl scaling the two theories, we fixed the expressions of 
the two scenarios in the Jordan and Einstein frames. This showed how they lead to markedly different phcnomenolo- 
gies, even if in the Jordan frame, they differ only by a constant term into the function multiplying the Ricci scalar: 
while the CQ case is bounded to produce effective corrections to gravity which increase its strength, the EQ case is 
able to provide corrections in both senses. This has important consequences which we highlight in the rest of the 
paper. 

We show that indeed linear perturbation growth may actually accelerate or slow down in EQ cases with respect to a 
case in which the dark energy is uncoupled, while in CQ scenarios perturbations can be enhanced only. In order to set 
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Correction type 


CQ 


EQ 


Gravity (AG/G.) DM-DM 


1 + 2C^ 


G/G. 


B-DM 


1 


G/G. 


B-B 


1 


G/G. 


Friction {AH/H) 


n 


1 


Mass (Am/m) DM 


^-Ca(4>-M 


1 


B 


1 


1 



TABLE L Summary of corrections required to run N-body simulations in CQ and EQ scenarios. 



up a suitable framework for the realization of N-body simulations in these cosmologies, we derive the Newtonian limit 
for both, again highlighting the different behavior that the gravitational corrections possess, never neglecting the role 
of Quintessence fluctuations. Finally, we list the quantities which need to be modified in N-body codes for taking into 
account the corrections to the expansion, strength of gravity as perceived by different species, as well as the modified 
perturbation growth rate, in body scenarios, when performing simulations of non-linear structure formation. 
Our work opens up the possibility of measuring the effect of non-minimal dark energy models in numerical simulations 
of structure formation, thus being able to constrain them through observations. N-body simulations nowadays reach 
the extension of thousands of Mpc, where indeed the interplay between large and small scale behavior of gravity really 
matters. Possible deviations from the ordinary behavior of gravity in General Relativity, which might be consistent 
with non-minimal dark energy models have to be detected coherently on all scales in order to provide convincing 
evidence. This capability can be built only by understanding properties and differences between the ways in which 
such deviations could manifest, developing capabilities for simulating and detecting them in large scale numerical 
simulations of structure formation and finally constraining them through observations. This work completes the first 
step in this path. 
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